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Abstract 



We consider CP violating effects in the decays r (Stt)//,- where both the 
J = resonance, ai, and = 0" resonance, vr', can contribute. The inter- 
ference between the ai and vr' resonances can lead to enhanced CP-violating 
asymmetries whose magnitudes depend crucially on the vr' decay constant, 
fjr'. We make an estimate of with a simplified chiral Lagrangian coupled 
to a massive pseudoscalar field, and we compare the estimates from the non- 
relativistic quark model and from the QCD sum rule with the estimate from 
the 'mock' meson model. We then estimate quantitatively the size of CP- 
violating effects in a multi-Higgs-doublet model and scalar-lepto quark models. 
We find that, while CP- violating effects in the scalar-leptoquark models may 
require more than 10^° r leptons, CP-violating effects from the multi-Higgs- 
doublet model can be seen at the 2a level with about 10^ r leptons using the 
chiral Lagrangian estimate of /^r' = (1 ~ 5) x lO^"' GeV. 
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1 Introduction 



The T has the same interaction structure as the e and the /x in the Standard Model 
(SM), apart from their masses. However, for practical purposes[0 the r lepton, the most 
massive of the known leptons, behaves quite differently from the e and /i leptons in that 
(i) the T has hadronic decay modesp], § (e.g. r ttu, pu, aiu, ...) which allow an efficient 
measurement of its polarization ^ Q and (ii) the couplings of the r to neutral and 
charged Higgs bosons |, 0] are expected to dominate those of the e and fi. These 
features make the r a rather special experimental probe of new physics. 

One phenomenon where new physics can play a crucial role is CP violation. Cur- 
rently CP violation has been detected only in the K meson system. The SM explains the 
effect adequately through a phase in the Kobayashi-Maskawa (KM) matrixp]. However, 
it is also possible that other sources of CP violation exist in nature. Recently the r decays 
into hadrons have been considered as probes of such a non-KM-type of CP violation in 
the scalar sector of physics beyond the SM. Nelson et a/.P] have considered the so-called 
stage-two spin-correlation functions to detect CP violation in the decay r {2tx)vt-, while 
in the same two-pion decay mode Tsai[0 has suggested tests of CP violation with lon- 



gitudinally polarized electron and positron beams at r-Charm factories. On the other 



hand, Kilan et a/.[ll| have studied T-odd triple momentum correlations in r decays into 
Kim and KnK. 

In the present work we consider the possibility of probing CP violation in the decay 
of the T into three charged pions, r (37r)z/^, where both the = l"*" resonance, ai, and 
= 0~ resonance, tt', can contribute. In particular, we investigate whether the large 
widths of these resonances in the decay of the r can be used to enhance CP-violation 
effects in extensions of the SM with scalar-fermion interactions which are consistent with 
the symmetries of the SM; the importance of broad resonances has been emphasized in 



the context of the top quark and the B meson [|T3[ in the last few years. 

In order to observe CP-violating effects there should exist not only a CP-violating 
phase but also processes interfering with different CP-conserving phases. In r decays 
one can in general have a CP-violating phase between the M^-exchange diagram and 
scalar-exchange diagrams in extensions of the SM such as a multi-Higgs-doublet models 
and scalar-leptoquark models. On the other hand, resonance enhancements of transition 
amplitudes and their coherent superposition can provide a large CP-conserving phase 
difference which leads to a significant enhancement of CP violating observables. 

The decay amplitudes of a r lepton into three pions, r {Stt)^^-, are dominated by 
contributions of the two overlapping resonances with different spins and relatively large 
width-to-mass ratios fl^ : 



ai : JP = 1+, m^, = 1.251 GeV, = 0.599 GeV, , . 

tt' : JP = 0", m^, = 1.300 GeV, T^, = 0.300 GeV. ^ ' 



Here we should note that even the parameters of the Oi are not accurately determined. 
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In the Svr decay mode of the r various phenomeno logical parametrizations p5|, |T6|, |17|, 
of the form factors have been used to analyze experimental datafl^, 0. We adopt the 
parametrization of the r-decay library, TAU0LA|1|], where the tt' parameters given in 
eq. (prrp can also be found. 

The three-charged-pion decay mode of the r is promising for the detection of CP 
violation due to the following reasons. First, no tagging of the other r is necessary. Second, 
the three-charged-pion mode can be measured not only at the conventional machines (e.g. 
CESR and LEP[0) but also at the planned B factories and at r-Charm factories where 
many r leptons (yearly 10^ to 10*) are expected to be produced. Third, reconstruction 
of the T rest frame is easy for the three-prong decay modes. Since at least two neutrinos 
escape detection it is in general difficult to reconstruct the r rest frame in r+r^ production. 
There are a few situations where the rest frame of the r can actually be reconstructed. 
One is r-pair production close to threshold where r leptons are produced at rest. This 
possibility can be realized at future r-Charm factories]^. Another is when both r leptons 
decay into hadrons. In this case impact-parameter methods allow us to reconstruct the 
rest frame of the r even for r's in flight. However, in the three-charged-pion decay mode, 
the direction of the r can be directly reconstructed through the precise and simultaneous 
determination of the r production and decay points. 

On the other hand, it is often difficult to make a reliable quantitative prediction 
for CP violation in hadronic r- decay modes due to the uncertainty in the hadronic matrix 
elements. This problem can, however, be removed by considering final states dominated 
by at least two neighboring resonances. The precise experimental determination of the 
widths and masses of the resonances will give a reliable handle on the calculation of 
CP asymmetries. In fact, the parametrization of the three-pion decay mode has been 
investigated quite thoroughly by many authors. Exhaustive discussion of the decay r — >■ 
SttUt- is given in Refs. p^ , p3| , 0. A complete parametrization of the ai/n' Stt decay 
currents can be found in the TAUOLA, which adopts from Ref. |]l6l /^/ = 0.02 ~ O.OSGeV 
However, the quoted value of /,r' may be overestimated because the mixing between the 
chiral pion field and a massive pseudoscalar qq bound state should be considered. Since 
any CP- violating effect in the decay r — > Sttz/^ depends on /^/, it is crucial to estimate its 
magnitude. We devote one section to this issue. 

The paper is organized as follows. In section 2 the decay distribution for r 
{3TT)i'r is presented. Its general form and our resonance parametrizations are given ex- 
plicitly. In section 3 CP-violating asymmetries in the decay r — >■ 3ni>r are introduced. 
In section 4 the multi-Higgs-doublet (MHD) model[^, Q and scalar-leptoquark (SLQ) 
models [^] are introduced as examples of models which can generate CP-violating effects 
in the decay r 37rz/^. In section 5 we estimate the magnitude of /t^' by employing a 
simplified chiral Lagrangian coupled to a massive pseudoscalar field. Finally, in section 
6 the possibility to detect CP violation in the three charged pion decay mode of the r is 
discussed quantitatively. Section 7 summarizes our findings. 
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2 Distributions 



The matrix element for the decay r (Svr) z/^ is written in the form: 
M = V2Gf [(1 + xMk, -)P+Yuip, a)J^ + rfu^k, -)PMp, ^)Jp\ , (2-1) 

with P± = (l±75)/2. Here Gp is the Fermi constant; p and k are the four- moment a of the 
r lepton and the tau neutrino, respectively, x ^-nd t] are complex numbers parametrizing 
the contribution from physics beyond the SM. The spin-quantization direction of the r~ 
is taken to be the direction opposite to the neutrino momentum (See Fig. 1(a).) and its 
helicity is denoted by o" (a = ±1). The tau neutrino is left-handed, so its helicity is — 1 
as indicated by u{k, —). J^^ and Jp are the vector and scalar hadronic matrix elements, 
respectively, and are given by 

J^=< (37r)- 1^7^(1 -75)u|0>, Jp =< (37r)-|rf(l + 75)M|0> . (2.2) 

Both the Oi and vr' resonances contribute to while only the pseudoscalar state, vr', 
contributes to J p. The explicit form of is found e.g. in the TAUOLA, and the hadronic 
(pseudo)-scalar current, Jp, is determined from the Dirac equation via the expression for 
J,: 

J, = N{^T^^\{q, - q,).F,{q\ s,) + {q, - q,).F,{q\ s,)] 

+q^CAsi{s2 - s^)F;{q\ s^) + s^i-Si - S;)F^{q\ s^)]], (2.3) 

Jp± = — - — q^Jf. 

vfiu + rud 



ml, 



rriu + rrid 



NC^,[s,{s2 - S3)F3(g^Si) + S2{s, - ss)F,{q^,S2)]. (2.4) 



Here qi and q2 are the four- momenta of two identical 7r~'s, q^ is the four-momentum of 
the 7T~^, and q is the four-momentum of the (37r)~ system in the decay of the r~. rriu and 
rrid are the current masses of the u and d quarks, respectively. The invariant mass squared 
of the 377 system, g^, and the three kinematic invariants, Sj [i = 1,2,3), are defined in 
terms of the three pion momenta, q^ {i = 1, 2, 3), as 

= (Qi + q2 + g3)^ -51 = (ga + g3)^ S2 = (gs + gl)^ ^3 = (gi + gs)^. (2.5) 

Also, T'^" = g^'^ — q^(f jc^ . The form factors, Fi {i = \ to 4), are normalized such that 
Fi(0, 0) = 1, and the coefficient 2v^/3 in ( p.3|) is fixed by the soft pion theorem. TAUOLA 
uses the parameters 

cos 6^(7 

= 7 ) ^n' = 7 2 ' y"^-^ 
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and adopts the parameter values 

cos^c = 0.973, 
U = 0.0933 GeV, 
dn'pTT = 5.8, gpj, 

It is convenient to cast the decay amphtude ( |2.1| ) into the form 



nip = 0.773 GeV, 

= 0.02 GeV, 
^ = 6.08. 



M = V2GFil + x) 



L A 



where the parameter ^ is given in terms of the x ^i-nd i] in eq. (|2.1| ) by 



v-x 



(m„ + md)mr \l + x , 



(2.7) 



(2.9) 



The r ^ Ut transition amphtudes, L^-a (A = 0, ±) and Lg-s, in eq. (|2.8|) are defined as 

= u{k, -)YP-u{P, (^KiQ, A), 
L,s=u{k,-)PMP,^), (2.10) 

where the polarization vector e(g. A) of the virtual vector meson ai satisfies 



(2.11) 



A=0,± 



The purely leptonic amplitudes are known functions of the kinematic variables and are 
given in the convention of Ref. |^ by 



0, 
rrir 



mt 



V2 



mt 



q^^o 



mt 



The hadronic matrix elements, H\ (A = 0, ±) and Hs, are 

2^2 



'-Ne,{q, A) [(g2 - qsTF.iq^ s,) + [q, - q,YF^{q\ s^)], 



Hs = NmrCn'[Sl{s2 - 53)^3(5^, Si) + S2(Si - S3)F4(g^ S2)]. 



(2.12) 

(2.13) 
(2.14) 



Assuming single meson dominance in each channel TAUOLA expresses the form 
factors Fi {i = 1 to 4) in terms of the meson propagators as 



Fiiq^s^] 
F3{q',s,] 



BAq^)Bp{s 



F^{q\s2) = BAq^)Bp{s2). 



(2.15) 



The ai, vr', and p meson propagators are parametrized ||14|| in the Breit-Wigner form with 
momentum-dependent widths: 



rrix — — imxTxiq"^) 



r 



X 



Jx{m] 



(2.16) 



for X = ai, tt' or p. The momentum-dependence of all widths has been determined from 
experimental data. The parametrizations of the widths available from the TAUOLA are: 



fAq' 



g2[1.623 + 10.38/g2 - 9.32/g4 + 0.65/g% for q^ > {nip + m^] 



4.1(g2 



(2.17) 



m 



9mly[l - 3.3{q - 9m2) + 5.8{q^ - 9ml 
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2\2l 



elsewhere, 



0, 

2 . 2 \ 3/2 
mi,— Ami. j ' 



0, 



for q > {nip + m^) 

elsewhere, 

for > (2m^)^ 

elsewhere. 



(2.18) 



:2.19) 



Ptt ill eq. ( |2.18| ) denotes the momentum of the tt in the (virtual) tt' rest frame. In addition 
to the parameters ( p..l| ) and ( |2.71) TAUOLA adopts for the p meson width 



= 0.145 GeV. 



^2.20) 



The amplitude for the r"*" decay into (37r) + , which is the CP-conjugated process of 
the T~ decay into (37r)~, can be determined in the same manner as that of the r~ decay 
into (37r)~. For the sake of discussion we use the same kinematic variables as in the r~ 
decay. In the r"*" decay qi and q2 are the four-momenta of two identical vr+'s, q^ is the 
four-momentum of the vr", and q is the four-momentum of the (37r)"'" system. The decay 
amplitude M for the process r+ {37!')^Ur can be written in the form 



M = V2Gf{1 + x) 



Y,L^xHx + {l + OL.sSs 



L A 



:2.2i) 



The r+ — * Ur transition amplitudes, L^x (A = 0, ±) and L„s, are given by 

= v{p, a)j''P.v{k, +)e;(g. A), 
Las = v{p,a)P+v{k,+), 



:2.22) 



where e(g. A) is the polarization vector of the virtual vector meson af with the helicity 
value A. We note that the hadronic amplitudes, Hx (A = 0, ±) and Hg, for the decay 
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r+ — s> (37r)+z/T- are the same as the hadronic amphtudes Hx{X = 0, ±) and Hs for the 
decay — (37r)^z/^: 

= ^^A, = H,- (2.23) 

Hence it is straightforward to obtain each amphtude of the r"*" decay into (37r)+ from the 
corresponding r~ decay amphtude through the CP relations among the amphtudes 
and the couphngs 

L^x = (-1)^L_,,_A, L^s = L.^s, (2.24) 
X = X\ V = r]\ (2.25) 

in the approximation where the imaginary parts of the intermediate W and scalar prop- 
agators are neglected. 

The decay of the t"^ into three charged pions depends on five independent kinematic 
variables. Because of two intermediate resonances, af and vr'^, decaying into (37r)^ it is 
convenient to consider both the rest frame and the (37r)^ rest frame. We define two 
coordinate systems as shown in Fig. 1. The two coordinate systems have a common y- 
axis which is chosen along the k x direction. Here k is the three-momentum of the 
tau neutrino. In the [x, y, z) coordinate system the z-axis is along the direction of the 
tau neutrino momentum, k, and the momentum is in the postive-x half-plane. The 
starred coordinate system is defined with a rotation by 6, the angle between k and ga, 
with respect to the common y-axis from the {x,y,z) coordinate system. We now consider 
the five variables; g^, si, S2, 0, and the azimuthal angle 0* of the vr^ momentum qi with 
respect to the x*-axis in the starred coordinate system (See Fig. 1(b).) 

The range of the variables is given by 

{3m^f < q^ < ml, {2m^f < < - m^f , 

S2min < S2 < S2max, < 61 < TT, <(/)*< 2tT , (2.26) 

where for given and si the minimum and maximum values of S2 values are 

1 



^2min ^ 

_ 1 
S2max 2 



+ 3m2 - si - ^(l-4m2/si) ■ \^l\q\ ml, s^) 



i + Zml - + V(l - 4m2/si) • \^l\q\ ml, si) 



;2.27) 



with the well-known triangle function, \{x, y, z) = x"^ + y"^ + z"^ — 2{xy + yz + zx). For 
simplicity we introduce the following notations: 



K = \{q\ml,Si) (^ = 1,2,3), G = V-A(Ai, A2, A3). (2.28 
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Using the kinematic variables and notations we obtain analytic expressions of the hadronic 
decay amplitudes (|2.14|) : 



N 



V2N 



Bai {'f') [A sin 6 + 5 (cos 6 cos 0* =F ^ sin ' 



Ba, (g^) [A cos9-B sin 9 cos 0*] 



(2.29) 

(2.30) 

(2.31) 



where 



A = SXslBpisi) + Bp{s2)] - (Ai - A2)[S,(si) - B,{s2)], 
B = G[B^{s,)-B^is2)], 

C = 2si{s2 - S3)Bp{si) + 2S2(S1 - S3)Bp{s2). 



(2.32) 



The hadronic amplitudes, Hx (A = 0, ±) and Hs, are symmetric under the interchange of 
two identical pions. The functions A and C also are symmetric, while the function B is 
antisymmetric. 

We denote the differential decay rates of the into three pions as 

d^ 

G{q^,Si,S2,cos9,(p*) = ———— ^j—T[t ^ (37r) z/^]. 



dg2dsids2d cos 9d(f)* 
d^ 

G(g^Si,S2,cos^,0*) = ^ ^77rr[r 

dg^dsids2d cos 



(37r)+z/^ 



(2.33) 



and examine the consequences of CP invariance in the differential decay rates. The 
CP transformation should relate the decay (37r)^z/T- to the decay r+ (371)"*" i/^. 

We find that CP invariance leads to the following relation in the differential decay rates: 



G{q^, Si, S2, cos6', 0*) ^ G{q'^, Si, S2, cos 6* 



(2.34) 



The relation enables us to construct a CP-conserving sum, S, and a CP- violating differ- 
ence, A, of the differential decay rates: 



S = G{q'^, si, S2, cos6', (f)*) + G{q'^, Si, S2, cos9, -c, 
A = G{q'^, Si, S2, cos9, 0*) — G'(g^, Si, S2, cos 9, — i 



(2.35) 



It is rather straightforward to obtain from the amplitudes (|2.8| ) and (|2.21|) the analytic 
forms of S and A; 



A 



1 + Re(0 Re{HoH: 



2F(g2)^[lm(0] [MHoH;) 



(2.36) 



where for notational convenience we have introduced the overall function F{q^) 



(2.37) 



From the expression ( |2.36| ) we note two important features of the CP- violating distribution 



A; (i) Every CP- violating asymmetry requires not only a non- vanishing Im(^) but also 
requires the interference between the longitudinal ai mode and the vr' mode, and (ii) the 
difference A depends only on cos0*, but it does not depend on sin0* which is odd under 
the naive time reversal T (See eqs. ( |2.3CI| ) and ( p.31| )). Here the transformation T means 
t — > — t without the interchange of initial and final states. 



3 Asymmetries 

As shown in section 2 every CP asymmetry requires not only a non-vanishing 
Im((^) but also the interference between the longitudinal ai mode and the pseudoscalar tt' 
mode; the interference is proportional to the /tt'. More explicitly A is proportional to the 
product of /tt' and Im(^), 

A cx [lm(0] ■ [/.'] . (3.1) 

In order to observe A it is useful to form an observable with an appropriate real 
weight function w(g^, Si, S2, cos 6*, 0*). A CP-violating scalar quantity is then obtained as 

{wA) = J i(7(g^, si, S2, cos 6', 0*) ■ A dg^dsids2dcos6'd0*, (3.2) 

where (X) means the integration of the quantity X over the allowed phase space of g^, 
Si, S2, cos6, and 0*. The statistical significance of this observable can be determined by 
the quantity 

- - (3.3) 



^(2) ■ (»'S) ' 

with S given in eq. (|2.36|) . The number of r leptons required to observe the effect at the 
1-cr level is then 

'V, = (3.4) 

where Br denotes the branching ratio of the r decay into three charged pions, which is 



6.8%|]T9l. By appropriately choosing the weight function, w, the CP asymmetry {wA) 
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can be made large. The optimal weight function maximizing the quantity e in eq. ( |3.3| ) 
is known p7| to be 



(2 a A*\ A(g^Si,S2,cos^,0 ) 

WopM , Sl, S2, cos^, (j) ) = — — . (3.5) 

S(g^si,S2,cos6',0*) 

We can also consider an observable with w = ±1 in phase space, which corresponds 
to the usual definition of an asymmetry. The distribution property of A in eq. (|2.36| ) sug- 
gests that we should consider two types of CP-violating forward-backward asymmetries, 
AipB and A2FB, whose weight functions are given by 

-Si, -52, cos 6', (p*) = sign[cos6'], 
W2(g^, si, S2, cos6', 0*) = sign[s2 - si] ■ sign[cos0*], (3.6) 

respectively. Note that the factor sign[s2 — si] is included in the definition of the A2FB 
asymmetry because the CP-violating asymmetry with only sign[cos0*] as a weight func- 
tion vanishes due to Bose symmetry. 



4 Models 

In the SM CP violation arises from a nontrivial phase in the KM fiavor-mixing 
matrix]^ in the hadronic charged current, but the KM- type CP violation can not be 
detected in r decays. As possible new sources of CP violation detectable in the r decay 
we consider new scalar-fermion interactions which preserve the symmetries of the SM. 
Then, it can be proven that only four types of scalar-exchange modelsPB| contribute to 



the decay r (37r)z/^. One of them is the multi-Higgs-doublet (MHD) model]^, 2^] and 
the other three models are scalar-leptoquark (SLQ) modelsl 



4.1 Multi-Higgs-doublet (MHD) models 

In this subsection we consider a MHD model with n Higgs doublets. The Yukawa 
interaction of the MHD model is 

Cmhd = QL^^^^dDR^ ++QlA'^^uUr^+Ll^F^^<^,Er^+ h.c. (4.1) 

Here Ql^ denotes left-handed quark doublets, and L^. denotes left-handed lepton doublets. 
Dr. (Ur.) and Er^ are for right-handed down (up) quark singlets and right-handed charged 
lepton singlets, respectively. The sub-index z is a generation index (i = 1,2,3). $j 
(j = 1 to n) are n Higgs doublets and $j = icr2$*. Sub-indices d, u and e denote the 
Higgs doublets that couple to down-type quarks, up-type quarks and charged leptons, 
respectively. and are general 3x3 Yukawa matrices of which one matrix can be 
taken to be real and diagonal. Since neutrinos are massless F^ can be chosen real and 
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diagonal. The MHD model has 2{n — 1) charged and 2n — 1 neutral physical scalars, 
and the Yukawa interactions of the 2(n — 1) physical charged scalars with fermion mass 
eigenstates read 



H+ + h.c. 



(4.2) 



Here Ma, Mu and Me denote the diagonal mass matrices of down- type quarks, up-type 
quarks and charged leptons, respectively. the positively charged Higgs particles. A^^^ 
are for left-handed neutrino fields and V for the KM matrix. Xj, Yi and Zi are complex 
coupling constants which arise from the mixing matrix for charged scalars. 

Within the framework of the MHD model, CP violation in charged scalar exchange 
can arise for more than two Higgs doublets ||30|, There are two mechanisms which give 
rise to CP violation in the scalar sector. In one mechanism [^, |3^ CP symmetry is main- 



tained at the Lagrangian level but broken through complex vacuum expectation values. 
However, this possibility has been shown to have some phenomenological difficulties |l33|, |6| . 
In the other mechanism CP is broken by complex Yukawa couplings and possibly by com- 
plex vacuum expectation values so that CP violation can arise from both charged scalar 
exchange and from exchange. CP violation in both mechanisms is commonly manifest 
in phases that appear in the combinations XY* , XZ* and YZ*. 

One crucial condition for CP violation in the MHD model is that not all the 
charged scalars should be degenerate. Then, without loss of generahty and for simplicity, 
we can assume that all but the lightest of the charged scalars effectively decouple from 
fermions. The couplings of the lightest charged scalar to fermions are described by a 
simple Lagrangian 



UiHD = n^Gp [x{UVMdDr) + Y{UrMuVDl) + Z{NlMeEr)\ H+ + h.c. (4.3) 
This Lagrangian gives the effective Lagrangian for the decay r SttUt 



£5f ^ = 2V2GeK 



X* Z Y* Z 

md^^{uLdR){VrLTR) + mu^^{uRdL){vTi^TR) 



M 



H 



H 



(4.4) 



at energies which are low compared to the mass of the charged Higgs boson. Then one 
can show that the contribution from the MHD model in the r Stti/,- decay distribution 



of eq. p.lj ) is represented by the parameters 

rrid 



Xmhd — 0, 



Vmhd 



Ml 



X*Z 



rridJ 



and CP violation in the MHD model is determined by the parameter 

rrid 



Im(e 



MHD) 



rriu + Hid. 



Ml, 



\m{XZ*) - ( ^ ) \m{YZ*) 
\mdJ 



(4.5) 



(4.6) 
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4.2 Scalar-leptoquark (SLQ) models 



In this subsection we discuss CP-violating effects from leptoquark exchange. There 
are three types of SLQ models which can contribute to the decay r dirur at 
the tree leveL The quantum numbers of the three leptoquarks under the gauge group 
SU(3)cxSU(2)lxU(1)y are 



$1 

$2 
$3 



3' 2. « 
o 



7 

6 
1 

3.1.-3 



3,3, — 
' ' 3 



[model I), 
(model II), 
(model III), 



(4.7) 



respectively. The hypercharge Y is defined to be Q = /a + Y. The Yukawa couplings of 
the leptoquarks to fermions are given by 



^SLQ = [-XijQL^En^ +x'ijUR^LL.]<!>i +h.c. 

^SLQ = h,QL,^r,Ll^ + y[pR.E'k^]^2 + h.c, 

4lQ = Zij[QL,lT2fLl] ■ $3 + h.c. 



(4i 



Here the coupling constants x^j , ypj and are complex when CP violation arises from 
the Yukawa interactions. Ql^ = {ui,di)L and Ll- = {i'i,ei)L- The superscript c denotes 
charge conjugation, i.e. ipR^L = "^T^T^V'S.l fo^^ ^ spinor field ip. f = {ti,T2,t^) and Tj 
{i = 1,2,3) are the Pauli matrices. In terms of the charge component of the leptoquarks, 
the Lagrangian relevant to the r — > SttUt decay is given by 



*-SLQ 

nil 
^SLQ 



f + h.c. 



4-^/'^ + h.c. 



(4.9) 



After Fierz rearrangement we obtain the effective SLQ Lagrangians for the decay r 



^11 



2Ml 



{dLUR){Ur^TR) + ^{dLCr^^URjiUr^af^uTR) 

{dLUR){f%u';^) + -^{dLa^''uR){f'j,a^,v';^) 



+ h.c.. 



1 2/13 1 



{dL-ff,UL){flj^KJ+ h.c 



Fl3| 



(47M«L)(^£7^^rJ+h.c. 



(4.10) 
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The tensor part as well as the scalar part in Model I and Model II can contribute to the 
r — s> (37r)z/^ decay. For simplicity, however, we consider only the scalar contribution in 
the present work. Then the size of new contributions from these three SLQ models is 



parametrized by the parameters (|2.1| ) 





0, 



















vLq = (4-11) 



2Mi 

2M'i 



^SL, = -1^, (4.12) 



VsLQ = 0. (4.13) 

Note that the vector-type interaction terms have only real coupling constants, and thus 
they cannot generate any CP violating effects. In particular. Model III does not contribute 
to CP violation in the decay r SttUt- so that the model will no longer be considered^. 
In Model I and Model II, the parameters governing CP violation are 

J ^ ml, Im[xi3x'i*3] 



(m„ + md)m, 4V2GfMI' ^ ' 

respectively, where all CP-conserving contributions from new physics are neglected in the 
normalization. This is justified because the contributions from new physics are expected 
to be small compared to those from the SM. 

4.3 Phenomenological constraints 



The constraints on the CP- violation parameters, ( [4. 61) and ( |4.14| ), depend upon 



the values chosen for the u and d quark masses. In the present work we use for the light 



u and d quark masses ||3^ 



m„ = S.OMeV, m^ = 9.0MeV, (4.15) 



and for the W boson mass we use Myy = 80 GeV. Inserting these values in (|4.6| ) and 
( [4.141 ) we obtain 



_^[lm{XZ*) - {5/9)lm{YZ* 



M^^hd) ^ -1.7 X 10-^^^ / )' \^ ' ^, (4.16) 



^However, in the semileptonic decays r KtitiVt and r KTrKur, the vector- type 
interaction terms from Model II and Model III in principle can have complex couplings 
to generate CP-violating effects. 
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Im(e| 



SLQ) 



SLQ) 



~ -7.5 X 10' 



My 

M, 



'Pi 



Im(xi3X 



My 

M, 



02 



Im(?/i3y;*3) 



(4.17) 



Clearly, sizable CP-violating effects require that lm.{XZ*) and ImiYZ*) are large and 
Mj^ is small compared to My/ in the MHD model, and similarly that \m.{xi^x'^^) and 
Im(|/i3|/5^^) large and M,^. (i = 1, 2) are small compared to Mw in the SLQ models. 

In the MHD model the strongest constraint on Im(XZ*) comes from the mea- 
surement of the branching ratio Br(i? — > XrVr) which actually gives a constraint on \XZ\. 
For Mh < 440 GeV, the bound on lm{XZ*) is given by 

lm{XZ*) < \XZ\ < 0.23M^GeV-^ (4.18) 

On the other hand, the bound P| on lm{YZ*) is mainly given by n^uu. The 

present bound is 

lm{YZ*) < \YZ\ < 110 for mt = 140GeV and Mh = 45GeV. (4.19) 

Combining the above bounds and assuming M^ = 45GeV, we obtain the bound on 
lm{^MHD) as 

IMa^D)! < 0.28. (4.20) 

The constraints on the leptoquark couplings can be obtained through several rare 
processes[^, ^ ^ As a typical bound we find from the experimental bound 

Bt[Kl —y fie) < 3.3 X 10~^^[^] for the branching ratio of the lepton-family-number- 
violating process — > fie, 



and from r(/iTi 



< 3 X lO-iiGeV-2, 

Ml 

eTi)/r(/iTi ^ capture) < 4.3 x 10"^^ 
\yuyl2 



Ml 



< 1.9 X lO^^^GeV 



r-2 



(4.21) 



(4.22) 



On the other hand, the helicity-suppressed vTej decay ^Q] can set a strong bound on 
and |yii|/n| under the assumption that Im(a;iix'{\) ~ Re(a;iix'{\) and lm{yiiy'{^) ~ 
^e{yiiy[\). Using the experimental value Re^p = (1.230 ± 0.004) x 10"^ [||] and the SM 
value R = 1.235 x 1O~^[0 for the ratio of B(7r — > ei^e(7)) to B(7r /iz/^(7)) we obtain 

|Im(xiia;;\)| \Re{xnXu] 



< 7 X lO^^^GeV 



Ml 
\lm{ynyu) 
Ml 



Ml 
Ml 



< 7 X 10-i^GeV-^ 



(4.23) 
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If we assume Im^xisx^^) ~ Im(x2ia;^2) ~ ^t^{xiix'*i) and lm{yi3y'{^) ~ Imd/ny^g) ~ 
lm{yiiy'{j) , then we obtain from ([4.21|) and ( [4.22|) the following constraints: 

\lm{esLQ)\ < 1.44 X 10-^ |Im(e|iQ)| < 0.90 x lO'l (4.24) 



Compared to the constraint ( [4. 201) on Im^^MHo), the constraints (|4.24|) on these 



SLQ CP-violation parameters are much more severe such that the SLQ models require 
a larger number of r leptons than are required by the MHD model to detect any effects 
of CP violation. However, we should mention here that the strong constraints ( |4.24| ) of 
the leptoquark couplings are based on the assumption that all the leptoquark Yukawa 
couplings are of the same size. This asuumption can fail in a class of leptoquark models 
where the lepton-family-number symmetry and the electron chirality are softly broken. In 
such leptoquark models one possible scenario is that the lepton-family-number-conserving 
couplings are much larger than the lepton-family-number-violating couplings, and also 
the couplings involving the third lepton generation (r and z/,-) are much larger than those 
involving the first lepton generation (e and z/g). In such a scenario the constraints ( |4.24| ) 
may be too stringent. 

One direct constraint on |xi3a;']*3| and l^/isyial "^^^ t)e provided through the mea- 
surement of B(r tiUt-). We assume that \\ii{xi^x'*-^) ~ Re(a;i3x']*3) and Imiyi^y'*.^) ~ 
Re(?/i3|/;*3). Then, employing the SM value P = (2.480 ± 0.025) x lO^^^ GeV[|T| and the 
experimental value Pexp = (2.605 ± 0.093) x 10"^'^ GeV||38|| for the r -nv^ decay width, 
we find that 

|Im(a:i3x;*3)| |Re(xi3x;*3)| n-6p.v-2 



|Im(t/i3Z/l3)| |Re(t/i3ya)| 



< 3 X 10-^GeV-l (4.25) 



It is clear that the direct constraints (|4.25| ) are by far weaker than the indirect constraints. 



( [4.21| ) and ( |4.22| ); therefore, the parameters Im(^5^Q) and Im(^|£Q) can in principle be 
significantly larger than the upper bounds ( [4.24|) . We expect the precise measurement of 
P(r tii^t) at future r-Charm factories to give a stronger direct constraint on |xi3a;']*3| 
and ly^y'ill- 



5 The n' decay constant 

As clearly shown in section 2 every CP-violating observable requires interference 
between the ai mode and the tt' mode, and thus the size of CP violation crucially depends 
on /tt' which determines the coupling of the tt' to the weak current. It is, therefore, 
important to make a reliable estimate of its magnitude. 

The most commonly used value of the vr' decay constant in the literature is obtained 
by a quark model calculation (the 'mock' meson model) by Isgur et al[^]. Especially, the 
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standard r decay Monte-Carlo program, TAUOLA[p!^, quotes from Ref. the value in 
the range 



= 0.02 ~ O.OSGeV, 



(5.1) 



where = 0.02 GeV is used as its default value. 

However, the decay constant of ir' is suppressed by the chiral symmetry which 
is not taken into account in the mock-meson calculation, and thus the value quoted in 
the TAUOLA might be overestimated. In this section we reconsider /^/ in the chiral 
Lagrangian framework and show that it actually vanishes in the chiral limit due to the 

/tt' is then proportional to the square of the 

HI 



mixing p2| of vr' with the chiral pion field, 
pion mass rnT^ and may hence be much smaller than 

In order to establish notations we start with the usual chiral Lagrangian with the 



flavor SU(2)lxSU(2)r symmetry[|3|; 



/: = /'trKxO + Ytr(x + X^), 

where the matrix-valued x ci-nd the Maurer-Cartan 1-form are defined as 



(5.2) 



with the chiral fields and given by 



t 



exp 



(5.3) 



(5.4) 



and tt" (a = 1,2,3) and T° are the Nambu-Goldstone fields and the flavor SU(2) group 
generators, respectively. The parameter b is given in terms of w-quark condensate (uu) 
and the vr decay constant / by 6 = —{uu)/p ~ 1.3GeV. The chiral fields and 
transform under SU(2)lxSU(2)p{ as 



(5.5) 



where the SU(2) group element h is introduced to maintain the relation i'ji = ^l- In the 



Lagrangian (|5.2| ) we have also introduced external scalar and pseudoscalar fields, s and p, 
and external left- and right-handed flavor gauge fields, and i?^^^. In the real world 
the external scalar and pseudoscalar fields are fixed to be s = diag(m„,md) and p = 0. 

We next describe how to introduce a massive pseudoscalar vr' in the chiral La- 
grangian formalism to discuss its physical properties. Since vr' is not a Nambu-Goldstone 
boson it should be treated as a matter field transforming as 



P hPhK 



(5.6) 
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where P is the vr' field. The low-energy effective Lagrangian of vr and vr' with the lowest 
derivatives can be written as 

Cp = /:chirai + ti{D^PD^P) - Mlti{P^) + iGpti{{x - X^)P) + 2Fptr(a^^D^P), (5.7) 

where Fp and Gp are parameters describing the low-energy properties of the tt'. The 
covariant derivative of the field P is given by 



D^P = d^P 



i\a 



Pi 



with the Maurer-Cartan 1-form 



a. 



{5.t 



(5.9) 



First, we discuss the chiral limit where x = 0- this limit the Lagrangian £p is 
expanded as 



Paf + fd^-K^A''^' - Fpd^P^A'^' + ■ ■ ■ , (5.10) 



where is an external axial vector field, A^^ = [R'^^ — Although the nonvanishing 

Fp in the last term of ( |5.1CI| ) seems to directly indicate the existence of the vr' decay 
constant the existence of Fp simultaneously causes a kinetic mixing of vr and P through 
the third term of ( p.lUD , and thus we need to diagonalize the kinetic mixing so as to 
evaluate the actual size of the vr' decay constant. 

We resolve this kinetic mixing by the redefinition of fields: 



\ 



(5.11) 



We find that the field redefinition pushes the vr' mass to a higher value. 



Mi 



2 ' 



and it forces the tt' decay constant to vanish in the chiral limit: 

U' = 0. 



(5.12) 



(5.13) 



Although our calculation is based on the effective Lagrangian with the lowest derivatives, 
the result ( p.l3| ) itself is rather general. Actually, we can show that the existence of the 
higher derivative terms do not change this result. 
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In the real world u and d masses are not zero and break the chiral symmetry 
explicitly. Including the vr-mass term due to the finite u and d quark masses and using 
the previous field redefinition ( 5.11 ) we obtain the mass term of the vr and vr' chiral 
Lagrangian as 

with ml = {rriu + md)b. We diagonalize the Lagrangian ( p.l4| ) in the perturbation of m^: 

Plugging (|5.15| ) into the Lagrangian ( ^.14| ) we find that the parameter S removing the 
mixing term is 



S-- 

The complete diagonalization leads to 



F| (Fp + 2Gp\ 



2 ^2 2 Ml, + AmlGpFp/p 

On the other hand, the meson coupling term to the axial vector current are written 
in terms of the physical fields TTpj^y^ and -Pphys as 



m. 



p 



and then the vr and vr' decay constants directly read 



/. = /, U = ^^/l-F'p/f{Fp + 2Gp). (5.19) 
As a result, we obtain the following expression: 

It is now clear that f^^/ is proportional to the square of the pion mass rrj-Tr, and thus it is 
naturally expected that the value of fj^' is small. 

The parameters Fp and Gp are expected to be of the size of QCD scale. We thus 
naively expect 

= 0(1). (5.21) 
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On the other hand, the value of the quoted in TAUOLA ( |5.1| ) leads to a large value 
of r^, 

= 300 ~ 6000, (5.22) 

which clearly contradicts our naive expectation (|5.21|) . 

Although the precise values of those parameters cannot be determined within the 
chiral Lagrangian framework, we can estimate them from the QCD dynamics. In the 
present work, we employ the non-relativistic chiral quark (NRCQ) model of Manohar and 
Georgi|^ and QCD sum rules|^5| to make a rough estimate of the parameters. 

In the NRCQ model the lowest S-wave bound states of the chiral quarks are iden- 
tified as p and it'. Neglecting the spin-dependent interaction we can take 

Mp ~ mp = 0.773GeV. (5.23) 

Both Gp and Fp are calculated from the wave function at the origin. In the non-relativistic 
approximation we find the relation 

2bGp ~ MpFp (5.24) 
On the other hand the value of Fp can be estimated from the relation 

< ^ 1 F2/,2 = (l•3GeV)^ (5.25) 

with = 0.0933 GeV. It is now straightforward to show that 

~ 4, ~ 2 X 10-^GeV. (5.26) 



Let us next discuss the evaluation of with the QCD sum rule technique ^5 



By using the operator product expansion of the isotriplet scalar current, and the 
pseudoscalar current, j^'^^ we find a sum rule, 

/ dse-'/^^ Imn(s) ~ 0, (5.27) 
Jo 

where n(s) is defined by 

n(-g2)5a^ ^ , I dxe^^^{T{ji-\x)jt\0) - ji^^^Hx)jl^'\0)}), (5.28) 

and M is the scale of the Borel transform. Assuming the Oq (isotriplet scalar meson), vr 
and tt' meson dominance we can derive an expression |[45|| for r^r; 

- - I l]exp[-{ml-ml,)/M% (5.29) 



— ml, 
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from ( p. 27] ) and its first derivative of Inserting the vr and vr' masses and mag = 0.98 

GeV, we find that 

~ 7, ~ 3 X 10-^GeV, (5.30) 

for the mass scale M around 1 GeV. 

Note that despite their very different approaches, the NRCQ model [^] and the 
QCD sum rule|^ give a similar estimate of the r^. According to those estimates one 
find f-j^i (2 ~ 3) X 10^'^ GeV, which is much smaller than the value (0.02 GeV) quoted 
in TAUOLA. Considering possible uncertainties in our estimates we use in our numerical 
analysis a slightly broader range of /,r': 

= (1 ~ 5) X lO-^GeV. (5.31) 

Certainly a more exact value of f-j^i can be determined directly from experiment. 



6 Numerical results 

For the numerical analysis we use for the CP-violation parameter values in the 
MHD model and the SLQ models: 



^HiMHo) = 0.28, (6.1) 
;|^q) = 1.44 X 10-^ Im(eliQ) 



Im(ekQ) = 1.44 X 10-^ Im(eliQ) = 0.90 x 10^^^ (6.2) 



which are maximally allowed values according to the constraints ( [4.20[ ) and ( [4.24[ ). We 



mention in passing once more that the constraints ( [4.24| ) on the leptoquark couplings may 



be too stringent in a class of leptoquark models where the leptoquark Yukawa couplings 
are significantly generation- dependent. For the tt' decay constant we consider the range: 
(1 ~ 5) X 10-3 GeV as in eq. ([CT| ). 



Fig. 2(a) shows the normalized differential decay width, E~^[dS/dv^g^, for the 
decay r (377)//,- as a function of the (37r) invariant mass, \fq^. Note that it has almost 
one peak. The reason is not only because m^/ and m^^ are quite similar in magnitude, 
but also because the tt' contribution is by far smaller than the Oi contribution. This 
justifies the approximation where all CP-conserving contributions from new scalar-fermion 
interactions are neglected in the normalization. In Fig. 2(b) we present the two CP- 
violating differential forward-backward asymmetries, dAipsld^f^ and dA2FB/dA/g^, as 
a function of the (37r) invariant mass In this case we take /tt' = 5 x 10"^ GeV and 
Im(^) = 0.28. The two asymmetries have a similar dependence on but it is clear 
that the asymmetry Aips is almost 4 times larger than the asymmetry A2fb- 

Table 1 shows the expected size of the integrated asymmetries Aifb and A2FB 
and the optimal asymmetry Sopt along with the number of r leptons, A^^, required to 
obtain the 2a signal in the MHD model and the two SLQ models for the CP-violation 
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parameter values ( |6.2| ) and for the vr' decay constant in the range ( |5.31| ). The optimal 
asymmetry, Sopu is optimized with respect to the kinematic variables (g^, Si, S2, cos6', 0*) 
so that its expected size denotes the maximally obtainable asymmetry in the three-pion 
decay mode. However, we find that the CP- violating forward-backward asymmetry Aips 
is sizable, while the other asymmetry, A2FB, is rather small. The Aifb asymmetry and 
the corresponding values of A^"^ given in Table 1 show that the CP- violating effects in the 
MHD model may be seen with less than 10*^ r's even for a rather small vr' decay constant. 
To see the CP-violating effects with the asymmetries AipB and A2FB in the SLQ models 
more than 10^^ r's may be required. On the other hand, the Eopt and the corresponding 
N'^ values show that the CP-violating effects from the MHD model can be detected 
with about 10^ r's, while still more than 10^° r leptons are required to see CP violation 
in the SLQ models. However, as previously mentioned, the CP-violating parameters 
lm{C,gi^Q) and Im(^|.^Q) may be considerably larger so that the required number of r 
leptons could be significantly reduced. In light of the fact that about 10^ and 10^ r leptons 
are produced yearly at B factories and r-Charm factories, respectively, we conclude that 
CP violating effects for the MHD model may be seen at the r-Charm factories through 
the CP- violating forward-backward asymmetry, Aifb, and even at the B factories with 
the optimal asymmetry, Eopt, for the n' decay constant in the range (|5.31|) . Finally, we 
note that, if we take /jr' = 0.02GeV, the value quoted in TAUOLA , the number of r 
leptons needed is reduced by a factor of at least 15. 



7 Conclusion 

In this paper, we have investigated the possibility of probing CP violation through 
the semileptonic decays r (Svr)//,- where two resonances, ai and vr', contribute. The 
interference of the ai and vr' resonances leads to enhanced CP-violating asymmetries 
whose magnitude crucially depends on the n' decay constant f^'- We made an estimate 
of f-rr' with a simplified chiral Lagrangian coupled to a massive pseudoscalar field, and we 
compared the estimates from the nonrelativistic chiral quark model and from the QCD 
sum rule with the estimate from the 'mock' meson model. We found that the mixing 
between the pion field and the massive pseudoscalar field renders /,r' proportional to the 
square of the pion mass, m^, and hence it is much smaller than the value quoted in the r 
library, TAUOLA. We considered two CP- violating forward-backward asymmetries, Aifb 
and A2FB, and the optimal asymmetry, Eopt, which is optimized with all the five kinematic 
variables (g^, Si, S2, cos6', 0*). With these asymmetries we quantitatively estimated the 
size of CP violating effects from the multi-Higgs-doublet model and the scalar-leptoquark 
models. We found that while the scalar-leptoquark models may require more than 10^° 
r leptons for the detection of any effect of CP violation, CP-violating effects from the 
multi-Higgs-doublet model may be seen at the 2a level with about 10^ r leptons using 
the chiral Lagrangian estimate of /^r' = (1 ~ 5) x lO^"' GeV. 
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Tables 



Table 1 The expected size of the CP-violating forward-backward asymmetries, AiPB {i ~ 
1,2), and the optimal asymmetry, Eopt^ and the number of r leptons, N'^ ^ needed for 
detection at the 2a level are determined for the range f^^i = (1 ~ 5) x 10^'^GeV with 
lm{iMHD) = 0.28 in the MHD model, and lm{^^sLQ) = 1-44 x lO'^ and lm(^|ig) = 
0.90 X 10-^ in the two SLQ models. 



Asymmetry 


Size(%) 








AiFB 


0.05 ~ 0.25 


12 ~ 300 


45 - 1125 


120 ~ 3000 


A2FB 


0.20 ~ 1.00 


0.75 ~ 19 


2.8 ~ 72 


7.4 ~ 185 


^opt 


0.48 - 2.4 


0.13 - 3.3 


0.49 - 12 


1.3 - 33 
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Figures 



Figure 1 The decay t~ — > {?>'k)~Vt viewed from the (a) t~ and (b) (37r)~ rest frames. The 
two frames arc related by a Lorentz boost along the tau neutrino direction and 
through a rotation by with respect to the common y-axis. 

Figure 2 (a) A plot of the normalized differential decay width, E^^fdS/dv^^, as a function 
of the (Stt) invariant mass, (b) CP-violating asymmetries as a function of 

the (Stt) invariant mass, \/^, with the tt' decay constant, f^' — ^ x 10~^GeV as a 
reference value. The solid line is for the CP- violating forward-backward asymmetry, 
dAiFB/d\/q^, and the long-dashed line is for the CP- violating forward-backward 
asymmetry, dA2FB/dy/q^. 
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